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In this work we describe, compile and generalize a set of tools that can be used to analyse 
the electronic properties (distribution of states, nature of states, . . . ) of one-dimensional 
disordered compositions of potentials. In particular, we derive an ensemble of universal 
functional equations which characterize the thermodynamic limit of all one-dimensional 
models and which only depend formally on the distributions that define the disorder. The 
equations are useful to obtain relevant quantities of the system such as density of states or 
localization length in the thermodynamic limit. 
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o ; Introduction 



The pioneering work of Anderson lH) changed completely the understanding of the properties of dis- 
ordered systems and meant the opening of a research field which is of primary importance nowadays. 
The physics of disordered systems is currently a significant part of condensed matter physics and it 
lO , has been the subject of an intense research activity specially during the last ten years. Electronic 
^ I localization due to disorder is a key element to understand different physical phenomena such as the 
I ' Quantum Hall Effect or the suppression of conductivity in amorphous matter. In the last years the ef- 
^ feet of the presence of statistical correlations in disordered systems has been analysed |l2l|3l|4l|5l|6l|7| 
^ and the conclusions regarding the appearance of extended states in the spectrum have been experi- 
^ , mentally confirmed for the case of short-range correlations [ 8 1 as well as for long-range correlations 
Q I ll9l[T0|. Scaling theory and Universality of the distributions of transport-related quantities characteriz- 
^ ' ing disordered systems are subjects which are still evolving nowadays: the conditions for the validity 
of single parameter scaling (SPS) have been recently reformulated [ ITl[l2|| and it has been found that 
different scaling regimes appear when disorder is correlated [ 13 1. The presence of disorder is of key 
importance for the characterization of low-dimensional structures, such as one-dimensional quantum 
Ch wires, since it plays a key role in the transport processes and it can strongly alter the electronic prop- 
^ ' erties of the system. Unlike the case of ordered matter, for disordered systems there is a lack of a 
general theory describing in a compact form their physical properties. Nevertheless a large ensemble 
>• . of different techniques exists that can be used to unravel some features of this kind of structures. 
^ ', Our fundamental premise to study the electronic properties of one-dimensional disordered systems 
I is to consider non-interacting spinless carriers within the independent particle approximation, that is 
- - ' the Hamiltonian of the system only includes the potential of a linear array of different atomic units 
. Also our approach focuses on the characterization of the static transport properties of these struc- 
tures. Within this framework, the aim of this work is to describe, compile and generalize a set of 
tools that can be used for all one-dimensional systems in order to analyse their electronic properties 
(distribution of states, nature of states, . . . ). In particular, we derive an ensemble of universal func- 
tional equations which characterize the thermodynamic limit of all one-dimensional models (within 
the approximations made above) and that are useful to obtain relevant quantities of the system such 
as density of states or localization length in that limit. Therefore a great part of our efforts are aimed 
at contributing to the growth of a general methodology that can be applied to all potential models in 
one-dimension. Let us also mention that the formalism here contained has already been used by the 
author and co-workers to describe successfully a large variety of one-dimensional disordered models 
Il35ll36ll37ll38ll39ll . However a complete and general derivation of the theoretical formalism is still 
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lacking; the present work comes to fill this gap. 

The work is organized as follows. In section [T] we make a thorough description of the continuous 
transmission matrix formalism and its applicability to finite-range as well as continuous potentials. 
Detailed analysis and calculations completing this section are contained in appendix|Xl The canonical 
equation and its derivation from the transmission matrix is treated in section |2l In section |3] the 
discrete transmission matrix formalism is briefly commented. The reliable parameters that can be 
used to characterize electronic localization are described in section IH where we particularly focus on 
the Lyapunov exponents, whose analysis is completed in appendix |Bl The procedure to calculate the 
distribution of states for the disordered chain is explained and generalize in section|5land appendix ICl 
to proceed subsequently with the construction of the functional equation formalism which is contained 
in section |6l and constitutes the main body of the work. The expressions of DOS and localization 
length in the thermodynamic limit in terms of the solutions of the functional equations are rigorously 
obtained. The applicability of the formalism developed is illustrated studying the one-dimensional 
tight-binding model. 



1 Continuous transmission matrix formalism 

The time-independent scattering process of a one-dimensional potential can be described using the 
well-known continuous transfer matrix method, 

A A fU,, M,A fAA _ ^ (AA ^ 



Br) \M2i M22; \Bl) \Bl^ 

where Al, Bl{Ar, Br), mean the amplitudes of the asymptotic travelling plane waves e'*^^, e~^^^, at 
the left (right) side of the potential. The peculiarities of the transmission matrix M and its elements 
depend on the nature of the potential. A detailed analysis on this subject can be found in appendix lAl 
As a summary let us say that for real potentials M belongs to the group SU{1, 1) and that the property 
det M = 1 holds for all kind of potentials whether they are real or complex. 
The transmission and reflection scattering amplitudes of the potential read 

1 r M21 R M12 

t = r = — r = — (2) 

M22' M22' M22' ^ ^ 

where the superscripts L, R, stand for left and right incidence. The insensitivity of the transmission 
amplitude to the incidence direction is a universal property. In general the reflection amplitudes will 
differ, although |r^| = |r^| for real potentials and complex ones with parity symmetry lfT4l . 

Obtaining the transmission matrix is specially easy for discontinuous short-range potentials such 
as deltas or square well/barriers, for which the asymptotic limit is not necessary to satisfy equation 
([T]). In these cases the effect of a composition of different potential units can be considered through 
the product of their transmission matrices, 

M = M7vM^_i---M2Mi, (3) 



therefore obtaining analytically or numerically the exact scattering probabilities of the whole struc- 
ture. This formalism can also be used to obtain the bound states from the poles of the complex 
transmission amplitude. An intuitive and general interpretation of the composition procedure can be 
given in the following form. Let us consider two finite range potentials Vi (x) , V2 (x) , characterized by 
the amplitudes ti, rf, rf , t2, r^, r^, and joined at a certain point. Then, the scattering amplitudes of 
the composite potential can be obtained by considering the coherent sum of all the multiple reflection 
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Figure 1 : A continuous potential 



processes that might occur at the connection region ifTSl . 




(4b) 



(4a) 



(4c) 



Replacing the scattering amplitudes with the elements of the corresponding transmission matrices 
Ml, M2, one can trivially check that in fact the latter formulae are the equations of the matrix product 
M2M1. Thus, the composition rules given by Q are not restricted to the convergence interval of the 
series ^^0(^2'^?)"- Th^Y provide an explicit relation of the global scattering amplitudes in terms of 
the individual former ones and can be easily used recurrently for numerical purposes. 

For continuous potentials the calculation of the transfer matrix is more complex. After solving 
the Schrodinger equation for positive energies, one has to take the limits x — » ±00 to recover the 
free particle states and identify the matrix elements. Hence equation ([T]) is strictly satisfied only 
asymptotically. However depending on the decay of the potential one could neglect its effects outside 
a certain length range. If the asymptotic transmission matrix of the potential in figured is 

known, then the matrix for the cut-off potential contained between the dashed lines can be written as 
(see appendix IXI) 



The cut-off matrix is the same as the asymptotic one plus an extra phase term in the diagonal elements 
that accounts for the total distance (^1 + ^2) during which the particle feels the effect of the potential, 
and also an extra phase term in the off-diagonal elements measuring the asymmetry of the cut-off 
((^2 — di). Doing such approximation one gets matrices suitable to be composed in linear arrays. 

2 The canonical equation 

As stated in the introduction we are treating one-dimensional atomic wires within the independent 
particle approximation. The electron-electron interaction is not considered and also the carriers are 
supposed to be spinless. Then, the Hamiltonian of the system only includes the potential of a linear 
array of different atomic units. From the solutions of the one-particle Schrodinger equation it is 
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Figure 2: Sites and sectors of a linear chain 
always possible to derive an expression with the following canonical form^ llT6l l4l 

vl>,+i = J(7,_i, 7.)^, - -^^vl/^-i, (6) 

where means the amplitude of the electronic state at the jth site of the wire, 7^ denotes the 
parameters of the potential at the jth site (jth sector) and the functions J(7j_i,7j), K{'jj), which 
depend on the potential and the energy, rule the spreading of the state from one site to its neighbours, 
as shown in figure El The canonical equation can be systematically obtained for a given solvable 
Hamiltonian and it contains the same information as the Schrodinger equation. It is not hard to see 
that J(7j_i, 7j) and K(-yj) can be chosen to be real functions provided the potential is real, so that 
the state amplitudes can also be considered to be real. Equation Q determines also the behaviour of 
other elementary excitations inside 1-D structures, thus it appears in different physical contexts such 
as the study of vibrational states (phonons), electron-hole pairs (excitons), . . . 

From the transmission matrix of the potential one can readily obtain the canonical equation apply- 
ing to the electronic states in the one-dimensional composite chain. Let us consider a linear composi- 
tion of potentials. All of them are formally described by the same transmission matrix with different 
parameters. And let be the transmission matrix of the jth potential. 



where the coordinates of the electronic wave function in the different sectors of the chain are chosen to 
satisfy that the amplitude of the state at all sites is simply given by the sum of the complex amplitudes 
of the travelling plane waves, that is = Aj +Bj for all j . To build the canonical equation one simply 
calculates the quantity '^j+i + X^j-i^ using and Mji^^ to write the amplitudes {Aj±i, Bj±i) in 
terms of {Aj,Bj). Then x is solved by imposing the coefficients of Aj and Bj to be the same. 
Following this procedure one concludes that the canonical equation for the most general potential can 
be written as, 

^,+1 = (s, + S,-i^^ - (8) 



where 



Sj = \ [(M,)ii + (M,)i2 + (M,)2i + (M,)22] , (9a) 
S3 = \ [(M,)ii - (M,)i2 - (M,)2i + (M,)22] , (9b) 
K, = i [(M,)ii - (M,-)22 + (M,)2i - (M,)i2] . (9c) 



^The meaning of the coefficients appearing in the canonical equation depends on the particular Hamiltonian. For a 
tight-binding model they have a straightforward interpretation in terms of the on-site energies and the transfer integrals, 
thus the equation is usually written in the form a^^'j = tj j+i^'j+i + For other models that comparison may 

not be so clear, so we keep a more general expression. 
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In the case of a real potential, using the symmetries of the transmission matrix (appendix |Al), one 
finds 



Sj = Re [(M,-)ii] + Re [(M,-)i2] , (10a) 
S, = Re [(M,)n] - Re [(M,)i2] , (10b) 
K, = Im [(M,)ii] - Im [(M,)i2] • (10c) 

And it also can be observed that for real and parity invariant potentials the functions Sj and Sj coin- 
cide because the off-diagonal elements of the matrix are pure imaginary. Then, the canonical equation 
can be easily calculated from the continuous transmission matrix of the compositional potentials of 
the system. 

Although the applicability of the canonical equation is not restricted by the ordering of the se- 
quence in the wire, it is a key ingredient to study non-periodic arrangements of potentials, for which 
the Bloch theorem is not valid. For certain boundary conditions, one can numerically obtain the per- 
mitted levels and the form of the envelope of the wave functions inside the system using equation 
Q. Apart from being useful from a numerical viewpoint, the canonical form also provides some 
analytical results concerning the gaps of the system's spectrum. For this purpose, the equation must 
be written as a two-dimensional mapping, originally proposed in reference ifTTI . that permits estab- 
lishing analogies between the quantum problem and classical dynamical systems lITSll . The matrix 
form of Q with the definitions Xj+i = ^j+i, Vj+i = "^j, reads 




(11) 



which in polar coordinates Xj = pj cos 6j, yj = pj sin 9j, leads to the following transmission relations 
for the phase and the moduli: 



-1 



= 7(9 j] 7j_i, 7j) = arctan <( ( J(7j-i, 7^) - tanOj ) }> , (12) 

^g^(^^.;^^_,,^^.)=cos2^,+ fj(7,_i,7,)cos^,-^^ . (13) 



Now let us impose hard-wall boundary conditions in our wire composed of atoms. That means 
^0 = "^N+i = 0. Using the mapping it is clear that the initial point is {xi,yi} = {\l'i,0} placed 
on the X axis. Thus for an eigenenergy, after all the steps the final point must be of the form 
{xN+i,yN+i} = {0, ^at} lying on the y axis. That means the whole transformation acts rotating 
the initial point. Therefore the permitted levels must be clearly contained in the ranges of energy 
for which the sequence of mappings generates a rotating trajectory (generally open) around the ori- 
gin, which is the only fixed point independently of the parameters of the mapping. This behaviour 
guarantees that after an arbitrary number of steps the final boundary condition could still be satisfied. 
However if all mappings have real eigenvalues the behaviour described is not possible (see for exam- 
ple reference llT9l ). And it follows that permitted levels cannot lie inside the energy ranges satisfying 



j2(7,_i, 7,) > ^Y(^y ^^"^^ 

Note that this conclusion does not depend upon the sequence of the chain, thus it holds for ordered 
and disordered structures. 
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3 Discrete transmission matrix formalism 



The problem of a one-dimensional quantum wire can also be treated via a composition procedure of 
another type of transfer matrices, when one obtains a discretized version of the Schrodinger equation. 
An analytical discretized form of this equation is given by the canonical expression Q, that can be 
written as 

The properties of the system can then be calculated from the product PatPat-i ■ ■ ■ Pi imposing ap- 
propriate boundary conditions. 

If the solutions of the differential equation are not known, one can always take a spatial discretiza- 
tion, translating the original equation 

il;"{x) = [V{x) -k^]ilj{x), (16) 

into 

V'n+l = {[Vn- k^] (AX)2 + 2}i;r^- (17) 

where we have defined tpn = i^in ■ Ax), Vn = V{n ■ Ax) and As being the spatial step. And the 
corresponding matrix representation is 

^n+A f[Vn - k'] (Ax)' + 2 -A / \ _ „ / \ „, 

Then, the scattering probabilities of the system can be numerically obtained by constructing Q = 
Q„Q„_i ■ ■ • Qi, considering a large enough distance n ■ Ax so that the correct asymptotic form of the 
state ■ip{x) = e''^^' + re^^''^ and ^(x) = te''^^ can be imposed at the extremes. The transmission and 
reflection probabilities are then given by 

^ 4sin^(fc- Ax) 

IQ21 - Q12 + Q22e'^-^- - Qiie-i'=-^-|' ' 



R 



hi - Q22 + Qi2e-^'=-^^ - Q2ie''=-^" 



hi - Q12 + Q22ei'=-^" - Qne-i'^-^- 
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4 Characterizing electronic localization 

The localized nature of the electronic states inside a disordered wire can be analyzed using different 
tools (see for example reference lEOll ). Let us see some reliable parameters which can be used as a 
probe of the localized or extended character of the carriers inside the system. 



4.1 Lyapunov exponents 

Lyapunov exponents emerge from random matrix theory lEDl . and they are used to characterize the 
asymptotic behaviour of systems determined by products of such matrices. They are a key element 
in chaotic dynamics [19| and play an important role in the study of disordered systems. For a full 
understanding of the meaning of the Lyapunov exponents and their expressions it is mandatory to 
recall Oseledet's multiplicative ergodic theorem (MET) (a complete analysis can be found in reference 
GJI), which in its deterministic version and without full mathematical rigour^ reads as follows. Let 
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be {M„} a sequence of d x d matrices and be Mn = MatMat-i ■ • ■ Mi. Then the following matrix 
exists as a limit ^ 

lim (M*^Mjv)^ = r ^ 0, (21) 

SO that its eigenvalues can be written as e'^^ < e'*'^ < ■ ■ ■ < e'*''' and the corresponding eigenspaces 
Ui, . . . ,Ud. And for every vector x of this d-dimensional space the following quantity exists as a limit 

A(x) = lim ^logllM^xll, (22) 

that verifies A(x) = max(Aj, . . . , A^) where {Ui, . . . , Uj} is the set of spaces in which x has a non- 
zero projection. The set {Aj} are the Lyapunov characteristic exponents (LCE) of the asymptotic 
product Mat. Therefore this theorem implies that the asymptotic exponential divergence of any 
spatial vector x under the action of the product of matrices Mat is determined by the LCE. More 
precisely, the divergence will be dominated by the component of x on {Ui} with the fastest growing 
rate. 

Now let us consider our one-dimensional quantum wires, which as already known can be de- 
scribed through products of different type of 2 x 2 matrices, namely the discrete transfer matrices Pj 
defined from the canonical equation in (fTSt and the continuous transmission matrices Mj defined in 
equation dl]). It can be proved that for one-dimensional Hamiltonian systems the two LCE come in a 
pair of the form {A, —A} (see appendix iBl. Considering the discrete transfer matrices we have 

xat+i = PatPjv-i ■ • -PiXi, (23) 
where xat+i ^ ( vl/ ) ^'^^ ^ I )' Therefore applying the MET 
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= aS^L N V + ^^^^ 

Imposing hard-wall boundary conditions A(xi) = limAr^oo ;^ log |\I/Af I which is straightforwardly 
equivalent to 

A= lim lyiog%:^, (25) 
j ■' 

a common expression found in the literature for the Lyapunov exponent, and that always provides the 
largest LCE 1 19|, in our case the positive one. 

On the other hand the same physical system can be realized using the continuous transmission 
matrix formalism, that must yield asymptotically the same values for the Lyapunov exponents if they 
have physical sense at all. Therefore, 

x^+i = M^M^_i ■ • • Mixi, (26) 
where xat+i = {^^] and xi = {^] corresponding to the amplitudes of the travelling plane 



waves. If we impose the initial conditions Ai = 1, Bi = r{E), then the final result will be A^+i = 
t{E), -Bat+i = 0, where r{E) and t{E) are the scattering amplitudes. Thus the MET implies 

-A= lim ^logr(E), (27) 



■'Several conditions must be satisfied by the set {M„} and its products that we suppose to be fulfilled in meaningful 
physical situations. 
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where T{E) is the transmission probability of the system, and obviously the negative Lyapunov ex- 
ponent is obtained. The above expression was first obtained by Kirkman and Pendry ll23l . It implies 
that for a given energy, the transmission of a one-dimensional disordered structure decreases asymp- 
totically exponentially with the length of the system Tn{E) ~ e^^'*'*^^)^ ll24ll . This is a consequence 
of the same asymptotic exponential decreasing behaviour exhibited by the electronic states for that 
energy. From this fact we define the localization length ^ {E) of the electronic state with energy E, if it 
exists, as the inverse of the rate of the asymptotic exponential decrease of the transmission amplitude 
with the length of the system for that energy, 

i{E) = \{E)-\ (28) 

This definition is also a measure of the spatial extension of the exponentially localized state inside 
the system, and it has a clear physical meaning. Although the Lyapunov exponent and therefore 
the localization length defined can only be strictly obtained asymptotically, it makes also sense to 
characterize the electronic localization in a long enough finite system through 

^(E)-' = -^\ogTiE), (29) 

because the Lyapunov exponent is a self-averaging quantity ll20ll . thus it agrees with the most probable 
value (its mean value) in the thermodynamic limit, for every energy. Therefore expression ( l29b gives 
relevant information of the localization length for finite N, since it will show a fluctuating behaviour 
around the asymptotic value. 

Finally let us say that a complex extension of the Lyapunov exponent is possible. 



Ac 

TV 



)-^El°g(^). (30a) 



-A,= lim ilog((B), (30b) 

N—i-oa iV 

\E'j being the complex amplitudes of the state and t{E) the complex transmission amplitude. The real 
part of this extension is related to the localization length whereas its imaginary part turns out to be tt 
times the integrated density of states n{E) per length unit of the system ll23l . 

4.2 Inverse participation ratio 

Alternatively, localization is also usually characterized by the inverse participation ratio (IPR) ll25ll . 
which is defined in terms of the amplitudes of the electronic state at the different sites of the system 
as 



V 1^ 



IPR = ' (31) 



For an extended state the IPR takes values of order N^'^ whereas for a state localized in the vicinity 
of only one site it goes to 1. The inverse of the IPR means the length of the portion of the system in 
which the amplitudes of the state differ appreciably from zero. 



5 Obtaining the density of states 

The density of states (DOS) g{E) gives the distribution of permitted energy levels and it is specially 
important for calculating some macroscopic properties of the structures which are usually obtained 
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from averages over the electronic spectrum. Strictly speaking g{E) is the function such that g{E)dE 
is the number of eigenvalues of the energy inside the interval (E, E + dE), and it is usually defined 
per length unit of the system. The integrated density of states n(E) is defined as 

n{E) = [ g{E')dE', (32) 



and measures the number of permitted energies below the value E. For a one-dimensional wire the 
integrated DOS can be determined from the imaginary part of the complex Lyapunov exponent. From 
(I30bl) one can write LI 6. .41 

n(E)^±^MnE)], „(^;) = J^,og(||j). (33) 

The electronic DOS can also be numerically determined using the negative eigenvalue theorem pro- 
posed by Dean for the phonon spectrum [26|, however this method cannot be applied for all potential 
models. It is possible to build a generalization of Dean's method to obtain the DOS for finite chains. 
This technique shows some relevant computational advantages comparing with the one involving 
the complex transmission amplitude. The whole derivation can be found in appendix O Defining 
Sjj+i = '^j+i/'^j the canonical equation ^ reads 

K('y) 1 

%,« = i(7,-.7,)-^— . (34) 

Now let us consider a wire composed of different atomic species {a} and let be y^a{E) the number 
of negative Sjj+i whenever 7^ = a, divided by the number of sites of the chain. That is, 'NaiE) is 
the concentration of a atoms after which the envelope of the electronic wave function with energy E 
changes its sign. Then the DOS per atom can be obtained as 



9iE) 



(35) 



Thus using the recursion relation (l34l) . one has to calculate the concentrations of changes of sing for 
the different atomic species, which must then be added or subtracted according to the sign of the 
functions K(a) for the corresponding energy. Finally a numerical differentiation with respect to the 
energy must be performed. 



6 The functional equation: the thermodynamic hmit 

The description of the properties of a system in the thermodynamic limit (TL) reveals the fundamental 
physics underlying the different problems, removing any accidental finite size effects. The TL tells 
us which observations are a consequence of a general physical principle. With this purpose Scaling 
Theory is intended to obtain different magnitudes in the TL by figuring out how they scale with the 
size of the system (see reference [27| for a thorough description of Scaling Theory) . Apart from 
Scaling Theory, a few authors have been in pursuit of obtaining analytically several quantities of 
an infinite one-dimensional disordered system. Dyson (1953) ll28l and Schmidt (1957)[29| derived 
analytically a type of functional equations for certain distribution functions containing information 
about the integrated density of states in the TL, for the phonon spectrum of a system of harmonic 
oscillators with random masses and the electronic spectrum of a delta potential model with random 
couplings, respectively. Although some efforts were made to solve numerically these equations 1.301 
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Figure 3: Sites and sectors of the one-dimensional system. 



|3Tl|32l|33l|34||, this approach was almost completely forgotten probably because of its cumbersome 
mathematics and the lack of analytical solutions. 

We assert that it is possible to derive a set of universal functional equations describing the TL 
of one-dimensional systems. In this way one can build a formalism which can be applied to a large 
variety of potential models. The solution of these equations can be used to obtain relevant magnitudes 
of the system such as the DOS or the localization length. Let us begin with the canonical equation 
describing our one-dimensional problem, 

= J(7,_i, 7.)^. - ^^^,-1, (36) 

where means the real amplitude of the electronic state at the jth site of the wire and 7^ denotes 
the set of parameters characterizing the potential at the jth site (jth sector). The functions J(7, P) 
and 1^(7) which depend on the potential and the energy, rule the spreading of the state from one 
site to its neighbours. From now on Greek letters are used to label the parameters of the different 
types of potentials composing the chain, while Latin letters always mean site indices. Using the map- 
ping technique described in section |2l one can define a phase 9 and radius p satisfying the following 
transmission relations: 

Oj+i = JiOj] 7j_i, 7j) = arctan | (^J (7j_i, 7^) - tan 6^^ | , (37) 

^) =:r(%;7,_i,7,)=cos2^^,+ (j(7,_i,7,)cos^,--^^ . (38) 

In order to ensure the continuity of the phase transmission for a given energy we work with the inverse 
function defined as 

7'^{6j+i + nn; 7j_i, 7^) = 'J"^(%+i; 7j_i, 7^) ± nn, Oj+i E [0, vr), n G Z, (39b) 

where the plus (minus) sign in ( I39bl ) must be taken when > (^-^^^^ < j corresponding 

to an increasing (decreasing) behaviour of the phase transmission. 

The goal is to calculate a distribution function for the phase 9, valid in the thermodynamic limit, 
so that the differential form of such a function acts as a natural measure of the phase in that limit. In 
this way one then would be able to obtain the thermodynamic average of any quantity of the system 
that could be written in terms of the phase. 
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The first step is to define the functions Wj{6) with 6 G [0,7r), that means the probability for 
9j ( mod vr) to be included in the interval [0,6'), that is dWj{9) means the probability that 9j ( mod tt) 
belongs to {9, 9+d9) , for a given energy. Therefore it follows that Wj (9) are monotonically increasing 
functions with 9 such that W^j(O) = 0, Wj{n) = 1 for all j. And we impose 



Wj{9 + nvr) = Wj{9) + n, 9 e [0,7t), n e Z, Vj. 



(40) 



According to the meaning of these distribution functions for the individual phases, it is clear that they 
must satisfy the relation 

dWj+^i9) = dW, (0-1(0; 7,_i,7i)) , (41) 

that is, the probability for 9j^i{modn) of being included in {9,9 + d9) must be the probability 
for 6'j( mod vr) of appearing in {7^^ (9), 7^^ (9) + d7^^{9)) ( modyr). Integrating the above equation 
leads us to 

Wj+ii9) = \W, (T-i(6^;7,-i,7,)) - W, (7-^(0; 7,_i, 7,)) | , (42) 

where the absolute value is necessary for the cases when 7^^ (9) decreases with 9 (i.e. [K{'jj) / K{^j^i) 
0), because the distribution functions must be positive. Since the inverse transmission function of the 
phase gives a value in the interval [—tx/2, 7r/2], the additional condition (l40b is used to ensure that 
the argument of Wj is always included in the interval [0, tt). And from the definition of the inverse 
transmission function it follows 



T-i(0;7,-i,7,) 



TT If -^(7j) 

2 X(7,_i 



if 



>0, 
< 0. 



(43) 



The equations relating the distributions for the phase at the different sites of the system clearly show 
that in fact those distributions only depend on the atomic species composing the chain. Thus the 
functions can be properly redefined in terms of the compositional species. Wj+i{9) is the distribution 
for the phase at the site (j + 1), generated after a potential of type 7-, (see figure |3l), therefore we 
relabel the function as Wj^i{9) = Wj^{9), that is the distribution function for the phase after a 7^ 
potential. And it is defined by 



where 



W,,{9)= W,^_,{7-U 
W^^i9 + n7r) = W^m+n, 



!7i-i,7i)) - W^7.-i (^) + ^{ij-mj] 
9 e [0,7r), n G Z, 



if 
if 



>0, 
< 0, 



(44a) 
(44b) 

(45) 



due to equations (l40b and (l43t . Now it is straightforward to carry out a thermodynamical average of 
the probabilities W^/. (9). We only have to sum over all the atomic species and binary clusters taking 
into account their respective concentrations. 



7 7:/3 



(46) 



where c-y is the concentration of the 7 species and C^p = Cpy is the frequency of appearance of the 
cluster -^P- or -/S'y-. Writing = c^p^p, where p^fs is the probability of finding a (3 atom besides 
a 7 atom, one can obtain an individual equation for each species, 

W,{9) = Y,p,p \Wp {7-\9-(3,^)) - Wp (I) + 5(A7) 

/3 

W^{9 + mx) = W^{9)+n, 



9 G [0,7r), n G Z. 



(47a) 
(47b) 



< 
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So that in the thermodynamic limit there exists a phase distribution function for each species com- 
posing the chain, and binary statistical correlations naturally appear in their definitions through the set 
of probabilities {pyp}. The completely uncorrelated situation corresponds to p^/j = C/3 for all species. 
Although we have supposed a discrete composition of the system, the same reasoning can be used for 
a continuous model in which the compositional parameters belong to a certain interval with a given 
probability distribution. 

Therefore solving equations (l47t . one would be able to calculate the average in the thermodynamic 
limit of any quantity of the system that can be written in terms of the phase 9 as long as it is a periodic 
function with period tt. The latter expressions are the most general functional equations valid for all 
one-dimensional systems for which a canonical equation of the form (l36b can be obtained. 



6.1 Calculating the localization length and the DOS in the thermodynamic 
limit 



Let us consider the Lyapunov exponent given by 

A 



(48) 



where ( ■ ■ ■ ) denotes the average in the thermodynamic limit and the amplitudes of the state at the 
different sites are considered to be real. Using the two-dimensional mapping defined previously. 



A 



+ lo; 



cos 



cos 6i 



log 



cos 9j+i 


cos 9j 


cos 9j+i 


cos9j 



(49) 



The middle term vanishes because the cosine is a bounded function that does not diverge as the length 
of the system grows. On the other hand the argument of the logarithm in the last term takes only the 
values ±1. Since log(l) = and log(— 1) = ivr it readily follows 



^log 



Re (A) : 
Im (A) = — i ( log 



Pj+i 
Pi 



(50) 
(51) 



From equation (l38b the average of the real part can be easily written using the distribution functions 
for the phase and therefore obtaining the inverse of the localization length ^(E), 



cos 9j^i 


cos 9j 


cos 9j^i 


cos 9j 



(52) 



which integrated by parts can also be written as 



e(i?)-^ = i5^c,p,^log3^(7r;7,/5)-^5^c,p,;3 / 



(53) 



On the other hand the imaginary part of the Lyapunov exponent increases by ivr every time the 
wave function changes sign from one site to the next one. Therefore by averaging equation (BTT) over 
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all possible species at the site (j + 1) when the jth species is a 7 atom, and dividing by vr, one obtains 
the fraction of 7 atoms after which the state changes its sign. 



- ( log 

TX \ 



COS Oj^i 


cos 6j 


cos Oj^i 


cos 6j 



J J / j+1 

/3 



TC 



dW^i9)\og 



cosT(6'; 7, P) \cos9\ 
_|cosT(^;7,/?)|cos^ 

■k/2 



(54) 



since the transmission function always returns a value in the interval [— 7r/2, n/2], where the cosine 
is positive. Thus it follows that [1 — lV^(7r/2)] is the concentration of changes of sign for the 7 
species, 'N-y^E) as denoted in section|5l And from equation (l35t the density of states per atom reads 



g(E) 



E 

7 



sgn[ir(7)]c. 



dW, (I) 



dE 



(55) 



6.2 Particular case: The canonical equation reads ^'j+i = 7(7^)^'^ — ^'j-i 

The functional equations (l47t can be considerably simplified depending on the particular model of 
the one-dimensional system. Here we consider one of the simplest forms for the canonical equation, 
appearing for example in the diagonal tight-binding model or the delta potential model with substitu- 
tional disorder. In this case the function 7(7) depends only on the parameters of one potential and one 
can take 1^(7) = 1 for all species. Therefore the problems concerning the changes of sign of the latter 
function are completely avoided and the inverse transmission function for the phase is an increasing 
function for all energies which depends only on one atomic species, (9; 7). Then equations (l47t 
read 



Wy{9 + mx) = Wy{9) + n, G [0, tt), n G Z. 



+ 1, 



(56a) 
(56b) 



If one further considers the case of uncorrelated disorder, that is p-^p = for all 7, P, then a global 
distribution function for the phase can be defined W{9) = J2-y c^Wy{9) being the solution of 



W{9) = J2c,W {7-\9;^)) - (|) + 1, 



W{9 



W{9)+n, 9e[0,n),nEZ. 



(57a) 
(57b) 



In this particular case only one functional equation needs to be solved and the localization length as 
well as the density of states per atom can be calculated respectively from 



aEy' ^ Re (A(E)) = ^ 5^ c, dW{9) log m 7), 

dl^(l) 



9iE) 



dE 



(58) 
(59) 
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7 Examples 



To exemplify the study of a quantum wire with the tools described in the previous sections, let us 
consider a basic one-dimensional model: a tight-binding Hamiltonian with nearest neighbour interac- 
tions, 

H = Y,{^k\^){^\+^k,k+i\k){k + l\+h^k~i\k){k-l\), (60) 

k 

where Sk are the energies of the on-site orbitals |A;) and tjj±i mean the transfer integrals, which we 
take equal to 1 for the sake of simplicity. The on-site energies follow a random sequence so that this 
model is said to have diagonal disorder. The one-dimensional Anderson model consists in choosing Ej 
from a finite continuous interval with a constant probability distribution. In our case the composition 
includes different discrete species {^i, £2, • • •} appearing with concentrations {ci, C2, . . .}. 

Since the orbitals { constitute an orthonormal basis of the Hilbert space of the system, the 
eigenstates can be written as |^) = Xlj^i The Schrodinger equation is then translated into a 
discrete equation for the coefficients Uj, 

Uj+i = {E - ej)uj - Uj^i, (61) 

showing the desired canonical form of equation ^ with J{sj) = E — £j and K{ej) = 1 for all j. 
Using the results of section |2l concerning the gaps of the spectrum, it is straightforward to conclude 
that the energy values satisfying \E — Sj] > 2 for all j, are not permitted. Therefore the eigenvalues 
can only be located inside intervals 4 units of energy wide centered at the different on-site energies. 
In fact each interval corresponds to the allowed band of the pure chain of each species. The simplest 
mixed system is a binary chain composed of two different species. In this case the spectrum only 
depends upon the quantity l^i — 62] thus the on-site energies are usually defined as ei = —e, 62 = £■ 
When 8^2 the eigenenergies are all included in the interval [—2 — e,e + 2]. That is the reason 
why this model is commonly referred to as a one-band model. If e > 2 a gap appears in the range 

[2-8,8-2]. 

From the canonical equation, the two-dimensional mapping defined in section |2l is easily built 
yielding the transmission functions 



1 

tan 6* . 

\2 „„„2 , 



T-i(0; Sj) = arctan ( E - 8j - ) , (62) 



3^(^; 8j) = 1-{E- 8j) sm{2e) + {E - 8^)' cos' (63) 

which only depend on one species at each step. This latter property together with the fact that K{8j) = 
1 for this model, simplifies considerably the functional equations (l47t . as described previously . For 
a chain with uncorrelated disorder a unique distribution function for the phase can be defined W{9) 
being the solution of 

W{9) = J2c^W {7-\9;8^)) - (^) + 1, (64a) 

7 

W{e + n7r) = W{e)+n, 9 E [0,7r), n E Z. (64b) 

Thus only one functional equation needs to be solved. And the DOS per atom as well as the localiza- 
tion length can be obtained in the thermodynamic limit from 

^E)-' ^ X{E) = IY.^, log 3^(7r; e,)-^^. ^7 £ ^(^^^^^^^ (65) 



9iE) 



dE 



(66) 
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Figure 4: log \t{E) \ vs length for disordered binary chains, for different energy values, on-site energies 
and concentrations. Only one realization of the disorder has been considered for each length. The 
coloured straight lines are plotted from the value of the Lyapunov exponent in the thermodynamic 
limit for each case. 




Figure 5: Comparison of the DOS for a fi- 
nite binary sequence with the limiting distri- 
bution for the infinite chain with parameters 
e = 0.5, ci = 0.3. 
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Figure 6: Lyapunov exponent for a binary 
chain with £ = 1 for different concentrations. 








Figure 7: DOS for a binary chain with e = 1 for different concentrations. The distributions for 
ci > 0.5 are the same but inverted with respect to the energy due to the symmetry of the Hamiltonian. 
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Using the discrete transmission matrix formalism, the scattering amplitudes of a finite chain can 
be obtained. In figure |4| the logarithm of the modulus of the transmission amplitude is plotted as 
a function of the length of the system for different binary chains. The exponential decrease of the 
transmission is clearly observed and the data for finite chains agrees with the value of the Lyapunov 
exponent obtained from the functional equation in the different cases. In figure |5l the DOS for a finite 
chain is plotted, showing a fluctuating behaviour around the distribution corresponding to the thermo- 
dynamic limit. The evolution of the DOS and the electronic localization length in the thermodynamic 
limit for a binary chain as a function of the concentrations can be seen in figures|6land|71 It can be ob- 
served how the tools used provide the correct results when the chain is composed of only one species: 
X{E) = and \{E) = arccoshdi? — e\/2) inside and outside the allowed energy band respectively 
and the density of states fits the correct form g{E) = vr"^ [4 — {E — e)"^] 

As can be seen in the analysis of the tight-binding Hamiltonian, taking the canonical equation 
as a starting point a systematic characterization of the electronic properties of the system can be 
performed, both for finite arrays via the transmission matrix formalism and in the thermodynamic 
limit with the functional equations. 

Aside from the tight-binding model, the functional equation formalism described in this work has 
been successfully applied by the author and co-workers to other one-dimensional potential models. It 
has been used to describe the DOS and the localization of the electronic states in the thermodynamic 
limit of a disordered array of delta potentials [35|, as well as to study the significance of binary 
correlations in the localization and transport properties inside the chain |36|||33. This formalism has 
revealed itself as a key tool to study other non-trivial 1-D models such as the Poschl-Teller potential 
13811391 which in the disordered configurations exhibits a large amount of exciting properties. We 
encourage the reader to check the given references since the results there contained are the best proof 
of the usefulness and power of this methodology. 

8 Final discussion 

In this work we have compiled and generalized some of the existing methodology to treat disordered 
quantum wires in one-dimension. The main result of the work is the derivation of a set of universal 
functional equations that describe analytically the thermodynamic limit of the disordered system inde- 
pendently of the potential model. The equations only depend formally upon the distributions defining 
the disorder in the system. The derivation of the functional equations has been described in a very 
detailed way, starting from the transmission matrix of the model and the canonical equation applying 
to the electronic states in the system. The functional equations can be solved numerically to carry out 
a systematic and thorough analysis of the electronic properties (DOS, electronic localization, effects 
of short-range correlations, . . . ) of several models of quantum wires in the thermodynamic limit, as 
it has already been done by the authors for different one-dimensional potentials Il35ll36ll37ll38ll39ll . 
The formalism developed in this work has made it possible to observe several interesting features of 
the electronic properties of one-dimensional disordered systems such as the fractal nature of the dis- 
tribution of states in the thermodynamic limit, the appearance of extended states or the drastic change 
induced by binary short-range correlations in the DOS and the localization of the electrons within 
the system. The functional equation formalism is not restricted to the framework of the electronic 
properties of disordered systems, since it can be applied also to other Hamiltonian systems, such as 
for example one-dimensional spin chains or one-dimensional phonons or excitons, which can be de- 
scribed in terms of a canonical equation of the form where the role of the electronic state is played 
by a different quantity. 

At the present time the theory of disordered systems does not include a unifying mathematical 
principle comparable to the Bloch theorem for the case of periodic structures. We would like to believe 
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that the derivation of the universal functional equations for disordered systems in one dimension may 
mean a little advance in this direction. The fact that independently of the potential model it is possible 
to write a general set of functional equations which only depend on the distributions defining the 
disorder and that characterize the thermodynamic limit of the systems, is to our minds a result that 
must be taken into account. And physical relevant quantities such as the DOS or the localization 
length can be directly calculated in the thermodynamic limit from the functional equations. In spite 
of their formidable aspect the equations may have analytical solutions in certain cases or they may be 
useful to extract analytically information of the systems in the thermodynamic limit. Work along this 
line will probably require the use of a tough mathematical formalism but it could also be very fruitful. 
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A The transmission matrix 

A.l Properties and symmetries of the transmission matrix 

Let be V{x) a finite range potential appreciable only inside the region [—d, d], so that the wave func- 
tion can be written as 

A2u{x) + B2v{x), -d<x ^d, (67) 

^^gifc(x-d) ^ ^^^-ik{x-d)^ x> d. 

u{x),v{x) being the linearly independent elementary solutions for each k of the continuum spectrum 
of the potential. By applying the continuity conditions of the state and its derivative at x — ±d, it is 
possible to reach an expression of the form 
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relating the amplitudes of the free particle states on the right and left sides of the potential. M is the 
continuous transmission matrix of the potential and its elements read in a general form: 

v{d)u'{-d) + v{-d)u'{d) - u{d)v'{-d) - u{-d)v'{d) 



Mil 



2W 

. kMd)v(-d) - k^ui-d)vid) + u'{d)v'{-d) - u'(-d)v'(d) 
+ 1 ^ ^ , (69a) 

_v{d)u'{-d) - v{-d)u'{d) - u{d)v'{-d) + u{-d)v'{d) 
2VV 

eu{-d)v{d) - eu{d)v{-d) + u\d)v\-d) - u'{-d)v'{d) 
+ ^ mv ' ^^^^^ 

v{d)u'{-d) - v{~d)u'{d) - u{d)v'{-d) + u{-d)v'{d) 

2VV 

^ . eu{d)v{-d) - eu{-d)v{d) - u'{d)v'{-d) + u'{-d)v'{d) ^^^^^ 

2kW ' 

v{d)u'{-d) + v{-d)u'{d) - u{d)v'{-d) - u{-d)v\d) 

JVI22 — 

2W 

k^u{-d)v{d) - eu{d)v{-d) - u'{d)v'{-d) + u\-d)v\d) 

+^ 2m ' ^^^^^ 

where W = v{x)u'{x) — u{x)v'{x) is the Wronskian of the solutions and it must be independent of 
X. A straightforward calculation of the determinant of the transmission matrix leads to 

v[—d)u'[—d) — u[—d)v'[—d) 

Therefore det M = 1 for all kind of potentials, since no specific assumptions have been made regard- 
ing V{x). Let us remark that this property is not a consequence of the time reversal symmetry of the 
Hamiltonian as it is usually stated. Let us study now the symmetries of the elements of the matrix in 
different special cases. 

A.1.1 Real potential G M) 

If the potential is real it is always possible to find real linearly independent solutions u{x)^v{x) for 
each value of the energy k. And from equations (l69b the following relations hold, 

M22 = MJ,, M21 = Mt2. (71) 

Therefore in the case of a real potential the transmission matrix can be written as 

M=f" f.V \a\^-W = l. (72) 



It is easy to check that these matrices satisfy 



M(i MM*=fi °). (73) 



-ly \Q -1 

The latter equation together with det M = 1 define the group SU{1, 1). 

V{x) with parity symmetry: If the potential is real such that V{—x) = V{x), then it is possible to 
find elementary solutions with parity symmetry. Let us suppose u{x) to be even and v{x) to be odd, 
then 

u{—x) = u{x), v{—x) = —v{x), (74) 

u'{—x) = —u'{x), v'{—x) = v'{x). (75) 
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Figure 8: Representation of the scattering process. 



Using this symmetry in (l69b one finds 

M21 = -M12. 

And together with conditions (1711) it yields a matrix of the form 



M 



a ib 
-ib a* 



b e 



\a\'-b' = l. 



(76) 



(77) 



A.1.2 Complex potential (r(x) G C) 

If the potential is complex, it is not possible generally to build functions u{x),v(x) being real, there- 
fore the conjugation relations (TtTT) are not satisfied. There are no special symmetries among the matrix 
elements. 



M 



a (3 
5 7 



07 — = 1. 



(78) 



V{x) with parity symmetry: In the case of a complex potential with parity symmetry, the same anal- 
ysis as for a real potential can be applied, and equation is obtained, since it does not depend on 
the elementary solutions being real or complex but only on their symmetries. 



M 



a 13 



a^ + P^ = 1. 



(79) 



V{x) with PT-symmetry: Let us consider a complex local potential invariant under the joint action 
of parity and time-reversal operations V*{—x) = V{x) [i40J . Then it is possible to find u{x),v{x) 
satisfying 



u [—X) = U[X), 



V [—X) = —v[x), 



Using the above symmetries in (l69b one is led to 



M2 



122 - ivx^i, 

Thus the matrix can be written as 



12 



-M 



12, 



M 



21 



-M21. 



M 



a ib 
ic a* 



b,c e 



Ur + fee = 1. 



(80) 
(81) 



(82) 



(83) 



A.2 Scattering amplitudes 

The scattering amplitudes are directly calculated from the transmission matrix. FigurelHis a pictorial 
representation of equation (l68t . Considering left incidence then Ai = 1, Bi = r^, A2 = t^, B2 = 0. 
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V{x) {x\V\x') = V{x)6{x-x') 



Transmission Matrix 



Real, 



All 
Real 



Complex 



Complex, 
Complex, 



V,V 

v,vf 



det M = 1 

a ih 



—lb a* 



a (3 

5 7 

a (3 

-P 7 

ic a* 



Scattering 
Amplitudes 



fj-i L J Fi, 



Table 1: Symmetries of the transmission matrix and scattering amplitudes for different potentials. 
Greek letters mean complex elements of the matrix while Latin ones represent real coefficients. V 
and T denote the parity and time-reversal operators respectively. 



And it follows 



1 



M. 



21 



M22 M22 

In the case of right incidence Ai = 0, Bi = t^, A2 = r^, B2 = 1. And the equations yield 



t 



R 



22 



(84) 



(85) 
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The insensitivity of the complex transmission amplitude to the incidence direction is trivially proved. 
Using the properties of the transmission matrix for particular cases of the potential is easy to see that 
for parity invariant potentials (real or complex) = r^, for generic real potentials = e^'^r^ and 
for generic complex ones both amplitudes differ. 

In table lA.2l a summary of the symmetries of the transmission matrices and scattering amplitudes 
is given for the different type of potentials described. 



A.3 Transmission matrix for a continuous potential 

For the most general continuous potential, equation (l68b is only satisfied asymptotically, that is the 
amplitudes of the asymptotic states 



^(-00) = Aie''^^ + Eie-''^^ 
^(00) = Ase''^" + ^se"''^", 

can be related via the asymptotic transmission matrix A^, 



As 
B3 



M 



(86a) 
(86b) 



(87) 



The procedure to obtain this asymptotic matrix is the following. First is solving the Schrodinger 
equation for positive energies so that the more general state reads 



^!{x) = A2u{x) + B2v{x) 



(88) 
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in terms of the elementary solutions u{x), v{x). Now one needs to build the asymptotic form of the 
elementary solutions, 

u{±oo) = f/i±e'^" + f/2^e-'^", (89a) 
i;(±oo) = l/i^e^^^ + l/a^e-''^^ (89b) 

Therefore the more general asymptotic state becomes 

^(±00) = {A2U^ + 52^1^) e''^" + {A^Ui' + ^2^^) e-''^". (90) 

Equating the coefficients with those corresponding to the asymptotic forms (l86b yields 

Ai = A2U{ + 52Vf , A3 = A2U+ + ^2^1+, (91) 

= A^U^ + B^V^- , B-i = A2U+ + B2V+ . (92) 

Solving (A3, in terms of (Ai, Bi) one obtains for the elements of the asymptotic matrix 

2i^ /f/+v- - v;+f/2- v^U{-utv,-\ 

W \UtV2 - V+U2 v+u^ - utv^~) ' "^^"^ 
where W = vu' — v'u is the Wronskian of the solutions. 

A.3.1 Including a cut-off in the potential 

Let us suppose that due to the nature of the potential, it is only appreciable inside the region [—di, ^2] 
(figure [l])- Then the transfer matrix for the potential with the cut-off relates the amplitudes of the 
plane waves dXx = d2 and x = —di, which can be written from the asymptotic forms (l86t . yielding 
the relation 

^^^^''''Vm (^^^-''''\ (94) 
B^^-ikd, I - ivicut I ^^^ikd, I , (y4; 



which implies 
leading to 



e 



ikd2 g \ /gifcrfi g 



Meat = ( g ^-i,,, ( g e-ifcd, 1 , (95) 



Once the asymptotic transmission matrix is known, the cut-off matrix is straightforwardly built. And 
for a given potential it is usually easier to calculate the asymptotic matrix than to construct directly 
the cut-off version using the continuity conditions at the cut-off points. 

B The Lyapunov exponents 

Let us prove that for a one-dimensional Hamiltonian system, suitable to be described in terms of 
products of random matrices, the two Lyapunov characteristic exponents (LCE) come in a pair of the 
form {A, —A}. 

Let us consider our system within the discrete transmission matrix formalism. If the Hamiltonian 
of the system can be written in terms of a potential V{x) then it can be shown that the canonical 
equation takes always the form 

^.+1 = ^(7,-1, 7,)^, - (97) 
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where 7^ denotes the parameters of the potential at the jth site of the system and J, K are functions 
depending on these parameters and the energy. Hence the discrete transmission matrix reads 



Pj(7i-i,7i) 








1) 



(98) 



The above matrix becomes a symplectic transformation if _ft'(7j)/i^"(7j_i) = 1 but in general det Pj ^ 
1. The necessary and sufficient condition to ensure that the LCE of the asymptotic product Pa? = 
Ptv • • • Pi (N — > 00) are {A, —A} is that the determinant of Oseledet's matrix T equals 1. Let us 



remember that T = lim 



Dt 



1 

I 2N 



. Thus we only need to prove that limAr_^oo(det P 



N 



,1/iV 



which is trivial if the individual transfer matrices are symplectic, like for example in the tight-binding 
model with constant transfer integrals or for the delta potential model, but it may not be so obvious in 
the most general case. For the proof let us suppose that our system is composed of two different kind 
of potentials 1, 2, in a random sequence. Then the product matrix will contain four different types 
of matrices, namely P(l, 1), P(l, 2), P(2, 1), P(2, 2) according to the different pairs of potentials 
occurring in the sequence and participating in the canonical equation. Then, 



lim (det P^l 



lim [(det P(l, l))^"(det P(l, 2)f'%det P(2, l))^2Hdet P(2, 2))^^^] 



(det P(l, l))^"(det P(l, 2)f'^{det P(2, l)f{det P(2, 2)) 



\Cl2 



C21 / 



\C22 



(99) 



Nij being the number of times that P{i,j) appears, and Cij the frequencies of appearance of the pairs 
—ij— in the thermodynamic limit. From (l98t it readily follows 



lim (det Pjv 

N^oo 



I N 





C12 


\K{1)] 


C21 


[ml 




[k{2)\ 





(100) 



A simple reflection symmetry argument requires that C12 = C21 so the above expression equals 1. 
And finally the eigenvalues of T must be of the form e'^, e~'^. 

This proof can be straightforwardly extended for a case considering k different potentials or for a 
continuous model with parameters inside a certain range with a given probability distribution. 

This result about the Lyapunov exponent for a one-dimensional system must be naturally expected, 
by the fact that the same physical problem can be treated through the continuous transmission matrix 
formalism and it must lead to the same results. Therefore considering the latter matrices which have 
always determinant unity, it is obvious that det F = 1. 



C DOS from node counting 

James and Ginzbarg obtained the expression of the integrated density of states (IDOS) of a linear 
chain of potentials in terms of the changes of sign of the wave function inside the different sectors 
of the system [41 J. Their reasoning is the following. Let us consider a binary wire composed of two 
types of potentials Pi and ^2- And let be [xj, Xj+i] the jth sector of the chain including a Pi potential. 
Then the elementary solutions for positive energy in this cell fi{x), gi{x), can be chosen to satisfy 

/i(x,) = l, f[{x,) = 0, (101a) 
gi{xj) = 0, g'iix,) = l. (101b) 

Now let us consider a certain energy E for which gi{x) has pi nodes inside the given sector (the first 
one at xj). Then for this energy, possible states include ^/'(x) = gi{x) + /i/i(x) with (pi — 1) zeros 
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and tp{x) = gi{x) — with pi zeros in the cell (for sufficiently small fi). For low E, pi = 1 and 

as the energy grows the index pi increases by one whenever gi{xj+i) = 0. Thus the energy spectrum 
can be divided in intervals according to the value of the index pi, so that if E lies in the interval 
labelled with pi then the solution ■ip{x) will have pi or (pi — 1) nodes in every sector of type Pi. To 
determine whether the number of zeros is pi or (pi — 1), one has to check if tp{xj) and ^(xj+i) have 
the same signs (even number of nodes) or opposite signs (odd number of nodes). Thus for that energy 
the number of nodes in a /3i sector can be written as 

pi-^ + (-ir|, (102) 

where 2 = 1 if ipi^Xj) and ^(xj+i) have the same signs and z = —1 otherwise. The same reasoning 
can be used for the species ^2- And for an energy E with the labels pi for the species Pi and p2 for 
P2, the total number of nodes inside the mixed system and therefore the IDOS per atom reads 



n{E) =ci ^pi + 




2 y V 2 y (103) 

-{-iy^j4,{E)-{-irJ42{E), 

where ci, C2, are the concentrations of the species and 'Ni^E), 'N2iE), are the concentrations of 
changes of sign for each species, that is the number of cells containing a certain species in which 
the state changes sign (i.e. iJj{x) has opposite signs at the beginning and at the end of the sector) 
divided by the total number of potentials of the system. To obtain the density of states one needs to 
evaluate n(E) in the interval {E, E + dE) in which the indices pi, P2, can be considered to remain 
fixed. Therefore the only quantities that can vary in the differential interval are the concentrations of 
changes of sign. Hence, 



9iE) 



dE ' ' dE 



(104) 



And this expression is straightforwardly generalized for an arbitrary number of species. The main 
result is that to determine the DOS correctly, depending on the energy range one has to sum or 
subtract the changes of sign of the wave function at sectors corresponding to different species. 

Now let us see how one can know the indices pi, p2, in a practical way. The system is completely 
determined by the canonical equation 

= J(7,_i, 7,)^, - J^^.-i- (105) 

The functions J, K, can be obtained in terms of the elementary solutions of the Schrodinger equation 
in each sector of the chain. Making use of equations (l69t and Q the function Ki^jj) can be defined 
as 

where /^^ (x),gj^ (x) are the elementary solutions in the jth sector [xj, Xj+i] with a potential of type 
7j. Imposing the additional conditions (IIOII) it follows for the case 7^ = Pi that K{Pi) = ~gi{xj+i). 
That is, the function —K{Pi) takes the same values that the elementary solution gi{x), verifying 
equations (IIOII) . would reach at the end of every Pi sector. Thus, whenever gi{xj+i) as a function of 
the energy changes its sign, and therefore the index pi increases by 1, K{Pi) also registers a change 
of sign. What is more, it is not hard to see that in equation (11041) the terms (— 1)^\ (— 1)^^ can be 
directly identified with the signs oi K{Pi), K{P2). And finally, one can write 



3(E) 



dE ■ ° ^ ^' -'^ 



(107) 
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From a numerical viewpoint one must do the transmission of the state through the system using 
the functional equation and count the number of changes of sign from site to site for the different 
atomic species, to perform finally a numerical differentiation with respect to the energy and sum or 
subtract the different contributions of the species according to the sign of the function K{(3) for the 
energy considered. 
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